Protecting entanglement in superconducting qubits 
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When a two-qubit system is initially maximally-entangled, two independent decoherence channels, 
one per qubit, would greatly reduce the entanglement of the two-qubit system when it reaches its 
stationary state. We propose a method on how to minimize such a loss of entanglement in open 
quantum systems. We find that the quantum entanglement of general two-qubit systems with 
controllable parameters can be protected by tuning both the single-qubit parameters and the two- 
qubit coupling strengths. Indeed, the maximum fidelity -Fmax between the stationary entangled 
state, poo, and the maximally-entangled state, pm, can be about 2/3 « max{tr(poopm)} = ^max, 
corresponding to a maximum stationary concurrence, Cmax, of about 1/3 ~ C(poo) = Cmax. This 
is significant because the quantum entanglement of the two-qubit system can be protected, even 
for a long time. We apply our proposal to several types of two-qubit superconducting circuits, and 
show how the entanglement of these two-qubit circuits can be optimized by varying experimentally- 
controllable parameters. 




PACS numbers: 85.25.-j, 03.67.Lx, 03.67.Mn 



I. INTRODUCTION 



Quantum information processing using superconduct- 
ing qubits (see, e.g., P, |3, d, Q) has made remarkable 
advances in the past few years. One-qubit and two-qubit 
quantum circuits (see, e.g., H, S B H, Ei) have 
been realized experimentally in superconducting systems. 
One of the most important issues in quantum information 
processing is how to couple two qubits, which has been 
widely studied theoretically and experimentalh; in super- 
conducting quantum circuits (see, e.g., [a, ml, [7 
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To couple two qubits, there are two types of approaches: 
(1) direct coupling and (2) indirect coupling. Examples 
of qubit-qubit direct coupling include, e.g., capacitively- 
coupled charge qubits [5j, iSj or inductively-coupled flux 
qubits d, H, [13. Examples of indirect coupling in- 
clude, e.g., qubit-qubit couplin g vi a a quantum LC os- 
cillator or an inductance Tsl. Il6j|. a Josephson junc- 
tion or an extra superconducting qubit acting as a cou- 
pler a nanomechanical oscillator , or a one or 
three-dimensional transmission line resonator [21, 22, 23]. 
The main merit of indirect coupling is that any two 
qubits can be selectively coupled in a controllable way 
(see, e.g., [l^l). By tuning some control parameters, one 
can continuously adjust the coupling strengths between 
qubits, which can be further used to switch between one- 
qubit and two-qubit operations. 
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Although maximally- entangled states can be prepared 
either by directly- or by indirectly-coupled two-qubit 
superconducting quantum circuits [25|, another equiva- 
lently important issue is to protect these entangled states 
in an open environment. Some entanglement-protection 
strategies have been proposed. For instance, it has been 
recognized (26l . [27| that a common dissipative environ- 
ment, e.g., a heat bath, would be helpful to protect 
the entanglement of two-qubit systems. Specifically, a 
common dissipative environment would lead to a collec- 
tive decoherence channel which may induce a so-called 
decoherence-free subspace [11] to protect special two- 
qubit entangled states. 

Even though a single collective decoherence channel 
could be used, in principle, to protect two-qubit entan- 
gled states, it is very difficult to obtain such a channel 
in experimental solid-state systems (e.g., superconduct- 
ing quantum circuits) . Independent decoherence channels 
usually lead to disentanglement [29^ . This loss of entan- 
glement cannot be recovered by local operations and clas- 
sical communications, when the two qubits reach their 
stationary states. Even worse, a mixture of the "inde- 
pendent decoherence" and the "collective decoherence" 
channels may destroy the decoherence-free subspace and 
lead to a failure of the entanglement protection. Thus, it 
is very challenging to protect entanglement in the pres- 
ence of such a dissipative environment, with independent 
decoherence channels acting on each qubit. 

In this work, we propose a two-qubit entanglement- 
protection strategy. We find that quantum entanglement 
can be partially protected even when the two qubits inter- 
act independently with their own separate environments. 
Our proposed strategy requires that both the single- 
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qubit oscillating frequencies and the coupling strengths 
between qubits should be tunable. Therefore, our pro- 
posal could be applied to several typical superconduct- 
ing quantum circuits, and other tunable solid-state qubit 
systems. 

This paper is organized as follows: in Sec. |TT] we 
present our main results for general two-qubit systems. 
The entanglement protection for directly- and indirectly- 
coupled superconducting qubits via an inductive or a ca- 
pacitive coupler is presented in Sec. IIIII and Sec. IIVI In 
Sec. |Vl we study the entanglement protection in super- 
conducting qubit circuits interacting with controllable 
squeezed modes in cavities or resonators (e.g., circuit 
QED). Conclusions and discussions are given in Sec. IVII 



II. GENERAL RESULTS 

We consider two coupled qubits with a general Hamil- 
tonian 
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where the Planck constant h is assumed to be 1. Here, 
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^ ±iay \ and ax' , o-y" , <a^z' are the ladder and 

Pauli operators of the j-th qubit. The frequency of the 
j-th qubit is denoted by u^aj- The real coefficients [i\ 
(with phase Q\) and [ii (with phase d-i) correspond to 
the a'^ a^-^ and o'L^Vj^'' coupling strengths, which are 
assumed to be tunable parameters. 

The qubits also interact with uncontrollable degrees of 
freedom in the environment. If the two qubits interact 
independently with their own environments, then, under 
the Born-Markov approximation [soj . we can obtain the 
following master equation: 

2 2 
p = -t\HA, p] + ^ V{D\a^l\ + ^ 2V^V\a^p\p, (2) 

where the super-operator 'D[L]p is defined as: 
V[L]p^LpL^~U^Lp~^pL^L, 

and Fi, represent the relaxation and pure dephasing 
rates for each qubit, respectively. In order to simplify 
our discussions, it is assumed that the two qubits have 
the same relaxation and pure dephasing rates. 
Below, we will use the concurrence C{p): 
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C{p) = max{Ai - A2 - A3 - A4, 0}, 



(3) 



to quantify the quantum entanglement between the two 
qubits (see Ref. [311), where the A^'s are the square roots 



of the eigenvalues, in decreasing order, of the matrix 

M = p(aWaf)p*(aW4^)), 

and p* is the complex conjugate of the density matrix p. 
If the two qubits do not interact with each other, i.e., 
= /^2 = in Eq. the stationary state of the two 
qubits should be the ground state p^ — |00)(00| (see, 
e.g., [HI), where the superscript "u" refers to "uncon- 
trolled" qubit system. Since p^ = |00)(00| is a separable 
state, the entanglement between two qubits is completely 
lost even if they are initially prepared in a maximally- 
entangled state. However, the following discussions show 
that the entanglement can he partially protected by tuning 
the interaction strength pi and the single-qubit frequen- 



in Eq. Q). 



A. Strong interaction regime: pi ^ 10 ^{ujai +t^a2) 

Let us now study the regime where the coupling 
strengths pi and p2 are about one-order of magnitude 
smaller than the sum of the two single-qubit frequencies: 
51 = bJai + L0a2, and the phases 0i, 6*2 are both time- 
independent parameters. For such systems, we have the 
following results: 

The solution p{t) of Eq. ([2]) tends to a stationary state 



Poo = VPm + (1 - P)Ps 



(4) 



as a convex combination of a maximally- entangled state 

Pm'- 

Pm = ^(lOO)-f e'C'i-^'lll)) ((00|+e-*(«i-^)(ll 

exp[-z(0i - (/.)] \ 











exp[i(6'i - (/))] 
and a diagonal separable state p, 

Ps = 



\ 



/3 




where 



The subscript "m" is an abbreviation of "maximally- 
entangled" , and the subscript "s" refers to "separable" . 
|0), |1) are the two eigenstates of a single qubit. The 
parameters p, <f> can be expressed as: 



2T2/T1 + (f72 + MTl)/Mpl 
SFa 
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FIG. 1: Schematic diagram of the energy transitions between 
the four eigen-states: |00>, |11>, |10>, and |01>. The sohd 
arrows denote the decays caused by independent relaxation 
and dephasing channels. The dashed arrows represent the 
coherent superpositions caused by the interactions between 
qubits (with interaction strengths fii and /U2, respectively). 



r2 in Eq. (O is the dephasing rate that is defined as: 

r2 -ri/2 + r^. 

The concurrence C of the stationary state poo (here- 
after called the stationary concurrence) and the fidelity 
F between poo and the maximally-entangled state pm 
(hereafter called the stationary fidelity) are given by: 



max ■ 



F{poo) = tr(p„i/9oo) 
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(6) 



The derivation of Eq. ([6]) is given in Appendix |^ 

Our calculations show that the stationary concurrence 
C{poo) and fidelity F{poo) are not affected by the in- 
teraction strength p2- Indeed, p2 induces a coherent 
superposition of the two eigen-states |01) and |10) (see 
Fig. [T]). These two states always decay to the two-qubit 
ground state |00), when subject to independent relax- 
ation and dephasing channels. Thus, p2 does not af- 
fect the stationary state poc- However, the interaction 
strength pi induces a coherent superposition of the two 
eigen-states |00) and |11). Of course, |00) is already in 
the ground state, while the state |11) can be partially re- 
covered from 1 00) by the coherent superposition caused 
by pi. Therefore, the stationary concurrence C{poo) and 
fidelity F{poc) only depend on pi. 

Equation ^ shows that the maximum concurrence 
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FIG. 2: Maximum COIlClirrGIlC6 Cmax cLllci maximum fidelity 
fmax, versus the ratio ri/r2 of the relaxation rates, as given 
in Eq. (0. Note that Cmax 0.31 and Fmax 0.65 when 
Fi/Fa ^ 2. 



C'max and the maximum fidelity -Fmax, where 



Cn 



Fr 




(7) 



can be obtained when the parameters pi and satisfy 
Fi + 6ATl 



Ml 



2 + + F2 



(8) 



The maximum concurrence Cmax and fidelity i^max , given 
in Eq. ([7]), are plotted in Fig. [2l This clearly shows that 
the concurrence Cmax and fidelity Fmax in Eq. ([7]) in- 
crease when the ratio 



(9) 



r2 ri/2 + r^ 

increases, and the highest concurrence and fidelity 

%/5- 1 



C„ 
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0.31, 
0.65 



(10) 
(11) 



can be obtained when F2 — > Fi/2, i.e., F^ 0. 

The obtained maximum concurrence Cmax niay not be 
high enough to be used for demanding tasks in quantum 
information processing. However, there are two possible 
ways to increase the stationary entanglement. First, the 
ratio ri/F2 in Eq. ([9]) can approach the optimal value 2 
for certain systems. For example, if the superconducting 
charge or flux qubits are at their degenerate points [l], 2 
3], the optimal value 2 might be obtained. In this case, 
the optimal concurrence and fidelity can be obtained, as 
shown in Eqs. (dH]) and (fTT|) . 
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Second, Eq. ([7]) shows that the maximum fidehty i^max 
between the stationary state p^o and the maximaUy- 
entangled state pm is always larger than 0.5, which makes 
it possible to introduce additional entanglement purifica- 
tion process to increase the proportion of the maximally- 
entangled state Pm- Even if the concurrence and fidelity 
reach the optimal values, as in Eqs. (fTU)) and pT|) . we can, 
in principle, further increase the stationary entanglement 
by using purification strategies, e.g., as in Refs. [3^ [33|. 



Weak interaction regime: pi <^ 10 ^{uJai + 



The results obtained in subsection III Al cannot be effi- 
ciently apphed to the case when pi, p2 10~^^2 (e.g., 
when pi, p2 are two orders of magnitude smaller than 
51). In fact, the optimal condition ([8]) shows that we can 
obtain the maximum concurrence Cmax and fidelity i^max 
only when 



Ml 



£1 

8 



> 



V2flrr+rf + r2 8(y5 + i)r2 



n. 



Thus, if pi ^ 10~^ri, then, optimally, the ratio ri/r2 
should be very small. In this case, from Eq. ([5]), the ob- 
tained maximum stationary concurrence Cmax will be ex- 
tremely small. Alternatively, in order to avoid this prob- 
lem, a time-dependent interaction between qubits should 
be introduced: 

If the phase 61 in Eq. ([1]) could be tuned to be 

01 ^ nt + tflQ ^ {uJal + LL>a2)t + 00, 

the obtained long-time state poo will be a time-dependent 
state 



Poo{t) 



2p 



-Pm{t) 



MiEi 



2pi + V1T2 



Ps 



as a convex combination of a time-dependent maximally- 
entangled state 



Pm{t) 



1 



-i(Ot+0o-f ) \ 



and a time-independent diagonal separable state ps'- 



/ 1-3/3 



Ps = 



/3 



where 



'1- 



8r2 



MiEi 



El \2pi+T{£2 



The corresponding stationary concurrence and fidelity 
can now be expressed as: 



C(poo) 
F{Poo) 



= max 



Ml (El - Ml) 



2pl- 
Pi{Ti - pi 
Apl 



E2E1 
1 

2E1E2 ^2' 



For the same reason discussed in subsection III Al the 
stationary concurrence C(poo) and fidelity F{poo) are not 
affected by the interaction strength p2- The proof of 
Eq. (in]) is similar to the proof of Eq. From Eq. (fT2|l . 
the maximum stationary concurrence and fidelity 




can be obtained when 



Pi 



T1T2 



(13) 



(14) 



A higher stationary concurrence and fidelity can be ob- 
tained, as in the strong interaction regime, by increasing 
the ratio ri/r2. 

Below, we apply the above results to several supercon- 
ducting circuits and discuss how their parameters can be 
varied so that the stationary concurrence and fidelity can 
be maxima. 



III. DIRECT COUPLING BETWEEN 
SUPERCONDUCTING QUBITS 

A. Two capacitively-coupled charge qubits 

Let us first study the superconducting circuit shown in 
Fig. m where two single Cooper pair boxes (CPBs) are 
connected via a small capacitor @1 ■ The Hamiltonian 
of the total system can be 



j=i 



''90 > 



- Ej{<^xj) cosf^j] + 4Jniri2, 



(15) 

where is a phase operator denoting the phase drop 
across the j-th CPB; fij = —id/{d(pj), which represents 
the number of Cooper pairs on the island electrode, is the 
conjugate operator of (pj. The reduced charge number 
Ugj, in units of the Cooper pairs (2e), can be given by 
Ugj = —CgVgj /2e, whcrc the parameters Cg and Vgj are 
the gate capacitance and gate voltage of the j-ih CPB. 
The Josephson energy Ej{^xj) of the j-th dc SQUID is 



i?j(<i>,,)==2i?0cOS TT-^ 




FIG. 3; Schematic diagram of two capacitively coupled 
CPBs. 



where Ej represents the Josephson energy of a single 
Josephson junction [s^l; ^xj denotes the external flux 
piercing the SQUID loop of the j-th CPB; and $o is 
the flux quantum. The coupling constant J between two 
CPBs is 

-r ^ 



where Cj and Cm are the capacitance of a single Joseph- 
son junction and the coupling capacitance between two 
CPBs. Ec = /2{Cg + 2C°j) is the single-electron charg- 
ing energy of a single CPB. For simplicity, we assume that 
Ec and E'j are the same for the two CPBs. 



Near Ug 



0.5, which is called the charge degenerate 



point, the two energy levels of the j-th CPB correspond- 
ing to rij =0,1 are close to each other and far separated 
from other high-energy levels. In this case, a single CPB 
can be approximately considered as a two-level system. 
In the charge basis, the Hamiltonian Ha in Eq. (fT5|) can 
be written as 0: 



r(2) 



(16) 

where Ec{ngj) — 4£'(7(1 — 2ngj) and the Pauli operators 
are defined as: 



a« = |0),,(l| + |l>„(0|, 

|o),.(o|-|i)„(i|. 



Here, |0)j and \ are the charge states with the Cooper 
pair numbers rij = 0, 1, respectively. 

Rewriting Eq. using the eigenstates of the single- 
qubit Hamiltonian, we have 



with Lo 



a] 



Ejy^\ Ec-j^EcM, andi^j, 



Ej{^xj)- The new Pauli operators ai^^ and al'" are de- 
fined by the eigenstates \+)j and of the j-th qubit 



l+>..(-| + (+1, 

l+>..(+|-|->..(-| 



where 



sin0j|l)j, 



|->j ^ sin6ij|0>j +cos6ij|l)j. 

Here, 0j = [a.Tcta.n{-Ejj/Ecj)] /2. 

Now, the two-qubit state p{t) evolves following the 
master equation: 



2 



2 



Let us now assume that the two qubits are both in the 
charge degenerate point Ugj = 0.5 with Ecj = 0, so that 
the dephasing effects can be minimized. In this case, we 
have = 0, which means that r2 = ri/2 -|~ P^ = ri/2. 
Further, from Ecj = 0, we have: 



(1) (2)^ 



(1)^(2) 



In some experiments (e.g., |5i]), the coupling strength 
J is of the same order of 



2 



n = ^2E°j cos ( vr^ 



$0 



where J « 4 GHz, E'j w 10 GHz and the dccoherence 
rates [sl Pi, r2 are of the order of 10-100 MHz < J, Q. 
If we let $2,1 = <I>3;2 = $a; and substitute 



n = 2Ej{^^), Pi 



J 

4 



into Eq. ([8|), then, in the limit ^ Pi, the maximum 
concurrence Cmax ~ 0.31 and fidelity -Fmax ~ 0.65, as in 
Eqs. pO)) and pT|) . can be obtained, if only 



J 



or, equivalently, 



V5-M 



Ej cos TT— 
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B. Two inductively-coupled flux qubits 



A(<I>'^) can be approximately expressed as [l7|: 



Let us now consider a superconducting circuit, as 
shown in Fig. IH where two flux qubits are coupled 
through their mutual inductance. Here, we modify 
the design used in the experimental device in Ref. Q. 
Namely, the small junction of each three-junction flux 
qubit is replaced by a dc SQUID, from which we can ad- 
just the tunnelling amplitude between the left and right 
wells (see Fig. [5]) of each single flux qubit. 




FIG. 4: Schematic diagram of two flux qubits coupled via 
their mutual inductance. 



Near the flux degenerate point « *&o/2, with the 
external flux (f>j piercing the superconducting loop of the 
j-th qubit, each flux qubit behaves as a two- level system. 
The total Hamiltonian of the two-qubit system can be 
expressed as 0: 



1 - 



$0 



where Lj is the self-inductancc of the superconducting 
loop of the j-th flux qubit, and 



cos 



is the tunable critical current of the j-th dc-SQUID with 
Iq being the maximum critical current. The coupling 
strength J between the two flux qubits is 

J = Mlpilp2, 

with the mutual inductance M between the two flux 
qubits. 




o,/o„ <l/2 



't,l<t>„ =1/2 



o ./o„ >l/2 



FIG. 5; Schematic diagram of the double-well potential of 
the j-th flux qubit with the two lowest-energy states for <l>j < 
$0/2, $j = ^0/2, and > $0/2. 

Let us now assume that $j = <&o/2, then the two-qubit 
Hamiltonian in Eq. (I17|) can be further simplified to 



where 



d'^f = \L,){R,\ + \R,){L,l 
Lj){L,\~\R,){R,l 



and \Rj) are the two lowest energy-states in the left 
and right wells of the j-th flux qubit (see Fig. [5]) . The 
parameter e(<&j) denotes the energy difference between 
\Lj) and \Rj) which can be expressed as: 



2I„ 



2 " 



where Ipj is the circulating current in the loop of the j-th 
qubit. The tunnelling amplitude A(<i>p between the two 
wells is tunable by varying the magnetic flux <i>J piercing 
the j-th dc-SQUID. In the limit 

0< ' ' - 1< 1, 

$0 



where 



and 



/3 V 



V2 
2 

V2 
2 



{\L,)-\R,)}, 
{\L,) + \R,)}. 



The decoherence process can be described by the master 
equation ([2]) under the Born-Markov approximation. At 
the degenerate point ($j = <I>o/2; j = 1,2), we have 
= 0, which means that T2 = {Ti/2) + T^ = Fi/2. 
Let us assume that $^ = $2 = and substitute 



n = 2A($'=), ^1 



J 
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into Eq. dH), then, in the Hmit A($'=), J > Ti, the 
maximum concurrence Cmax ~ 0.31 and fideUty i^max ~ 
0.65 can be obtained when 



75 + 1 



In experiments [9|, A($'^) and J are of the same order 
(~ 1 GHz) that are far larger than the decoherence rates 
Ti, r2 ~ 1-10 MHz (see, e.g., [sl]). Thus, the above 
optimal conditions could be reahzed in experiments. 





IV. TUNABLE COUPLING BETWEEN 
SUPERCONDUCTING QUBITS: 
STRONG-INTERACTION REGIME 



FIG. 6: Schematic diagram of two charge qubits coupled via 
an LC oscillator. 



There are two ways to tune the system parameters to 
achieve the optimal condition ([5]). One way is by tun- 
ing the sum of the single-qubit oscillating frequencies Q, 
which was used in Sec. IIIII In this section, we study an- 
other way to achieve the optimal condition ([5]) by tuning 
the coupling strength jii between the two qubits. 



A. Variable-coupling between two charge qubits 

Many strategies have been proposed to obtain a con- 
trollable coupling between qubits (see, e.g., p^. ITgI. [TtI. 
m, m [13, [B, [H, m, HI). Let us first study the su- 
perconducting circuit shown in Fig. [6l where two CPBs 
are coupled via an LC oscillator. This strategy was first 
proposed in Ref. [l3| and also investigated by other re- 
searchers (e.g., in Ref. 



V - o , |22|)- In the charge deg enerate 

point, the two-qubit Hamiltonian in Refs. [l|, [lal is 



(18) 



with Eint = Ej{<^^i)Ej{<^,2)IEl and Ej{<^,,) = 
2^;° cos(7r$a;j/$o)- Here 



|0), 



-*|0)„-(l|+z|l)„( 



and |0)j, |l)j are the two charge states of the j-th CPB. 
The quantity E^ in the expression above for the coupling 
strength E^^t can be written as: 



Er 



C. 



qb 



0_ 



where C^b = 2CjCg [2Cj + Cg] ^ is the capacitance of 
a single CPB in the external circuit, and L is the induc- 
tance of the coupling current-biased inductor. 



Rewriting Eq. (llSp under the eigenstates of the single- 
qubit Hamiltonian, we have: 



Ha = 



where 



Eint , (1) (2) , (1) (2)n 



a ^ 



(i) _ 



^ X 1 



and 



l+). = ^{|o),-|i).}, 

h>. = ^{|0). + |l),}. 

Now, let <^,j-i = $3,2 = ^x- Since r2 = ri/2 at the charge 
degenerate point, then, replacing and /ii in Eq. ([8]) by 
2i?j($a;) and i?i„t/4, the maximum concurrence Cmax ~ 
0.31 and fidelity F^ax ~ 0.65 can be obtained when 



cos TT 



-if 



1 



Et 



V5 + 1 EY 



f / V u r i^j 
when these conditions hold: E^^ E'j ^ Fi, F2. 



B. Variable-coupling between two flux qubits 

Our strategy can also be applied to flux qubits with 
controllable coupling. Here, let us consider a supercon- 
ducting circuit design as in Ref. [18] (see Fig. [7]), where 
two three-junction flux qubits (qubits 1 and 2) are cou- 
pled via an auxiliary three-junction flux qubit (qubit 3). 
This middle flux qubit (qubit 3), acting as a coupler, is 
connected to qubits 1 and 2 by sharing junctions a and b 
with the same Josephson energy Ej, while junction c is 



8 



smaller than a and 6, with Josephson energy ai^j, a < 1. 
By adiabatically eliminating the degrees of freedom of the 
auxiliary flux qubit 3, the total Hamiltonian of the flux 
qubits 1 and 2 becomes [l^ 



1 2 



(2) 
z ' 



J = l 



where e($j), Aj, j = 1,2, ai-''' and ai"''' have the same 
meaning as in subsection IIIII B. When a ^ 1, the cou- 
pling strength J($3) between the two flux qubits 1 and 
2 becomes [18|: 

From Eq. (|19p . the coupling strength J($3) is tunable by 
varying the flux $3 piercing the superconducting loop of 
the auxiliary qubit 3. 



^(3) 




FIG. 7: Schematic diagram of two three-junction flux qubits 
coupled via an auxiliary flux qubit. 

At the flux degenerate point, i.e., $j = <&o/2 for 
both qubits, and using the eigenstates of the single-qubit 
Hamiltonian, we have: 



1 2 



0) 



2 



^ J(|3) (^(1)^(2) +^a)^(2)). 

Since the two flux qubits are at their flux degenerate 
points, then, replacing f2 and /ii in Eq. ([8]) by Ai -|- 
A2 and J($3)/4, the maximum stationary concurrence 
Cmax ~ 0.31 and fidelity i^max ~ 0-65 can be obtained 
when J($3) « (Ai + A2)/(v^ + 1), i.e.. 



cos 27r 



$1 



^ (A1+A2) 

>oy " (V5-i-i)Jo' 

when these conditions hold: Ai, A2, Jo 3> ri,r2. 



From the experiment [T^, where Ai, A2 and Jo are 
of the same order (~ 1 GHz) and far larger than Fi, F2 
(~ 1-10 MHz), the above optimal condition could be 
satisfied by varying the magnetic flux $3 through the 
middle superconducting loop. 



V. TUNABLE COUPLING BETWEEN 
SUPERCONDUCTING QUBITS: WEAK 
INTERACTION REGIME 

In this section, we study how to protect quantum en- 
tanglement in superconducting circuits where two charge 
qubits are coupled to a one-dimensional transmission line 
resonator. Since the interaction strength (10-100 MHz) 
between the two charge qubits coupled via the resonator 
is far smaller than the single-qubit oscillating frequency 
(5-15 GHz), then we are now considering the weak- 
interaction regime. From the analysis in subsection III Bl 
in order to protect entanglement in this case, the follow- 
ing time-dependent interaction Hamiltonian should be in- 
troduced: 

The main idea in this section is the following: borrow- 
ing strategies that produce controllable squeezed fields 
in optical cavities (see, e.g., [38]), an auxiliary flux qubit 
circuit is introduced to squeeze the oscillating mode in 
the resonator (see Fig. . The auxiliary flux qubit cir- 
cuit in fact acts like a A-shaped three-level atom which 
is further driven by a classical field. By adiabatically 
eliminating the degrees of freedom of the auxiliary flux 
qubit circuit, one can obtain a controllable squeezed field 
in the resonator where the squeezed coefficient is tunable 
by changing the coupling strength between the classical 
driving field and the auxiliary flux qubit. With the help 
of the controllable squeezed field in the resonator, one 
can continuously adjust the stationary entanglement be- 
tween the two qubits. 



Controllable squeezed electric field in a 
transmission line resonator 



We first show how to obtain a controllable squeezed 
electric field [sol . [iol [4ll . [4^ in the resonator by using a 
theoretical proposal of realizing squeezed states in cavi- 
ties (ssj . The auxiliary flux qubit circuit in our proposal 
(shown in Fig. [S]) acts as a three- level system with A-type 
transition [43,, ^] . 

As depicted in Fig. [9l we are now considering a three- 
level system with a ground energy level \g), an interme- 
diate energy level and an excited energy level |e). 
Here, the transitions \g) ^ \i) and |e) ^ \i) are cou- 
pled dispersively to the quantized cavity mode in the 
resonator, with coupling strengths Xg and Ag. The tran- 
sition \g) |e) is coupled dispersively to a classical 
field with coupling strength Ad and frequency Q. In the 
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auxiliary flux qubit 



charge qubit 1 



X i X 



VgiO 



X 



X ®(t >< 



charge qubit 2 



X i X 



one can obtain the following reduced Hamiltonian by adi- 
abatically eliminating the degrees of freedom of the three- 
level artificial atom |38j : 

= Lu.a^a + ^(e-'(^*+'^»)a^2 ^ g^(m+0o)^2^^ ^21) 

where 



is the effective frequency of the cavity mode; ^ and (f)o 
are t/ie effective amplitude and the initial phase of the 
squeezed field. The relation between ^ and 4>o is given by 



transmission line resonator 

FIG. 8: Schematic diagram of our proposal for protecting 
quantum entanglement in two charge qubits coupled to a res- 
onator. 



Notice that one can continuously adjust ^ by varying the 
coupling strength between the classical field and the 
three-level artificial atom. 



auxiliary flux qubit 

:0, 



X 



®o X 




transmission line resonator 
(a) 



FIG. 9: Schematic diagrams for realizing a controllable 
squeezed electromagnetic field in a resonator coupled to an 
auxiliary fiux qubit. (a) The auxiliary flux qubit and the 
transmission line resonator: the parameters in the auxiliary 
flux qubit are the same as those in Ref. [31; by varying the 
flux "I>s threading through the SQUID loop one can obtain a 
A-shaped three-level artificial atom, (b) Transition energy- 
level-diagram of the A-shaped three-level artificial atom. 



rotating-wave approximation, the total Hamiltonian of 
the three-level artificial atom and the resonator can be 
expressed a,s H = Hq + V, with 

Ha = i:)ca^a-iljc\g){g\+5\i){i\+LUc\e){e\, 
V = (Aga|i)(g| -l-h.c.) -f (Aea|e)(i| +h.c.) 

+(Ad|e)(g|e-'"*+h.c.), (20) 

where h.c. means Hermitian conjugate; Qc is the fre- 
quency of the resonator; 5 is defined as a detuning from 
the energy levels |e) and \g) to the intermediate energy 
level \i). 

Let us initially prepare the artificial atom in the in- 
termediate level \i). With the help of the dispersive- 
detuning condition: 

5» lAJ, lAel, \n-2CjX 



B. Tunable coupling between qubits 

In Fig. [51 let us now consider the interaction between 
the two charge qubits and the cavity field. After elimi- 
nating the degrees of freedom of the auxiliary three-level 
system, we can obtain the following total Hamiltonian of 
the charge qubits and the cavity field: 



H 



-j(Ot+0o)gt2 



a') 



+ ^ g{rij — cos ttjcr^^^ 4- sin aja^''){a) -f a) 
1 ^ 

+ ^'^^aj<J^^'^ +i^ca'^a, (22) 



where 



has the same meaning with the corresponding quantity 
in Sec. [m] and [IVl g = -e(Cg/Cs)Kms is the coupling 
strength between the resonator and a single qubit; Cs 
is the total capacitance of a qubit; V^^^ = yJujj2(X- is 
the root mean square (rms) of the voltage across the LC 
circuit; Cr is the capacitance of the resonator; r/j = 1 — 
2ngj; and the angle aj = aT:ctan[Ej{^xj) / Ec{l — 2ngj)]. 

By using the rotating-wave approximation and assum- 
ing that Ugj — 1/2 {j = 1, 2), the total Hamiltonian H 
in Eq. ([22l) can be rewritten as [30l |: 



_|_^(g-i(nt+0o)Qt2 gi(Ot+0o)^2-j_ 

Here, to simplify our discussions, we have set 
= 'i'x- 
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We now assume that the qubits and the cavity field are 
in the dispersive regime, i.e., 

A = [Ej{^^) - uic] - 100 MHz > I5I - 10 MHz. 

Thus, we can introduce the foUowing unitary transfor- 
mation to diagonaUze the Hamiltonian Hjc- 



U = exp 



IB 



Up to first order in g/A, we have 

UHjcU^ « oj.a^a + ^e-*(^"'^*+^«)at2 ^ ^g,(nt+0o)^2 




2g^e 



-i{nt+<po) 



A A2 J 

where 



di^'^ + h.c. 



<:^a = £^j($.)+4gVA, 

By adiabaticaUy ehminating the degrees of freedom of 
the resonator, the following two-qubit Hamiltonian can 
be obtained: 



Ha 



■U) 



I (1) (2) 



E 



Here, we have omitted all the single-qubit terms induced 
by the interaction between qubits and the resonator be- 
cause of the conditions: 

i?,7(<i>,)/2»5VA, e^/A. 

As analyzed in subsection A of section|Vl we can continu- 
ously adjust the parameter ^, thus the coupling strength 
/ii is continuously tunable. 

Since the two superconducting charge qubits also inter- 
act with the uncontrollable degrees of freedom in the en- 
vironment (e.g., quantum noises induced by charge fluc- 
tuations on the electric gates), the discussed two-qubit 
system should be considered as an open quantum system. 
For this two-qubit system, the master equation ^ can 
be obtained under the Born-Markov approximation [2^ . 

From Eq. p^ . at the charge degenerate points for both 
qubits, we know that the optimal concurrence Cmax ~ 
0.31 and fidelity i^max ~ 0.65 can be obtained when 



1 



V5 + 1 g- 



-ri. 



Using now the same experimental parameters from 

Ref. pj: 

A = Ej -ujr^^ GHz - 4.8 GHz = 200 MHz, 
g = 20 MHz, ri/27r - 0.1 MHz, 
Ag, Ae - 10 MHz, 5 - 100 MHz, 

the squeezed amplitude ^ is of the order of 1 MHz, which 
can be realized by a strong microwave driving field with 
coupling strength A^ ~ 100 MHz. These parameters 
show that our entanglement-protection proposal is ex- 
perimentally realizable. 

Although this section mainly concentrates on how to 
protect entanglement in two charge qubits coupled to a 
transmission line resonator, our proposal is also extend- 
able to two flux qubits in a coplanar transmission line 
resonator [isj. Since the system parameters [45| are al- 
most of the same order of those for charge qubits, there is 
no essential difference between them, as far as applying 
our proposal. 



VI. CONCLUSIONS 

In summary, we propose a strategy to protect quan- 
tum entanglement against independent quantum noises 
in several types of superconducting circuits. For two 
superconducting qubits coupled strongly via an induc- 
tive or a capacitive element, one can tune the stationary 
entanglement by varying the single-qubit oscillating fre- 
quencies or the coupling strengths between the qubits. 
For superconducting qubits weakly coupled via a quan- 
tum cavity (e.g., a transmission line resonator), an auxil- 
iary superconducting three-level system [i^, [4J| with A- 
shaped transition is introduced to induce a controllable 
squeezed field in the cavity. Such a controllable quantum 
squeezed field can be further used to entangle two qubits 
in open environments. Optimally, one can obtain a max- 
imum concurrence of about 0.31 and a maximum fidelity 
of about 0.65 for the above systems. 

Even though the proposed strategy can be used to pro- 
tect quantum entanglement, the obtained entanglement 
may not be high enough to be used in quantum infor- 
mation processing. Additional entanglement purification 
processes (e.g., [SJl) should be introduced to increase 
the stationary entanglement. These procedures could 
make the superconducting circuit too complex. For this 
reason, further research (possibly using the methods in 
Ref. [401) will be focused on modifying our proposal to 
obtain higher stationary entanglement. 

Another interesting problem would be to develop 
a short-time regime entanglement-protection strategy, 
e.g., to protect entanglement during the gate operation 
process. In this regime, the correlation effects of the 
environmental noises should be considered, which leads 
to non-Markovian noises [i^l. Different effects may be 
produced under non-Markovian noises. The existing 
decoherence suppression strategies [i^, against non- 
Markovian noises may be helpful to solve this problem. 
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where $7^ (z = 1, • • • , 15) are all traceless basis matrices 
in (jAip and TOj = tr(f2ip). 

Let m = (mi, • • • , mi^)^ , and thus the master equa- 
tion (21) can be rewritten as [13, HH : 



(A2) 



rh = Oa m + D m + g, 



where Oa is the adjoint representation matrix of —iHa 
and (Dm + g) is the coherent vector representation of the 
Lindblad terms: 

2 2 

with D < and g a constant vector. Further, let 



APPENDIX A: DERIVATION OF THE 
MAXIMUM CONCURRENCE AND FIDELITY 



(m-^4, TO^4, TO23, TOI3) , 



to'' = (mi4, 77123, TO^^)^, 



In this appendix, we show the derivation of the con- 
currence C of the stationary state poo and the fidelity 
F between poo and the maximally- entangled state pm- 
Thus, we will derive Eq. ^ in the main text. 

In order to simplify our discussions, let us use the so- 
called coherent vector picture as in Refs. [13, [EH. Con- 
sidering the inner product (A, Y) — tr(A^y), we can find 
the following matrix basis for all two-qubit matrices: 

y iji) 1^(1) 



where 



{1 T oy c>y 1 /-r 

2-'4x4, ''I4i "l4; "23' ^'23' 2 ■ 



23' 2"^ ' 

1 (2) 1 (2) 1 (1) (2) 1 (1) (2) 1 (1) (2) 



ia«4^\l7f„17i3,^-i'^-i'^}' (Al) 



where /4x4 is the 4x4 identity matrix, and r2i4, ^\!^-, 



^23, ^23j ^14' ^23 are defined as: 



"14 



^'23 



x/2 



f7f 



Vo 
V 



1 

%/2 



1 

V2 



7f\ 



v4 / 



"14 



^23 



V5 \ 



V \/2 



V2 



Vo 

r 

V 



V2 



1 

%/2 



~1 
\/2 



0/ 



The system density matrix p can be expanded under this 
matrix basis as: 

1 

p = ^-^4x4 + 'y^miVtj, 



"^14 = tr(ri!j»4p), 771^3 = tr(f7^3p), a,/? = x,y, z, 



ruaf} = tr 



a,/3 = 0,a;,y,z, 



and CTq"'-' = /2x2 (i = 1, 2) are 2x2 identity matrices 
acting on the qubit j . Then, we can rewrite Eq. (jA2[) as: 



777' 



777 



1=1 

4 

J^M-O'^^W + D'^m'^ + g'\ (A3) 

i=l 

4 

^7t,0,'777' -Fi:»'777% 



i=l 



where L)p = -4(ri + 2r^)J4x4 = -8r2/4x4 and 

7*1 = cos 01, 7i2 = 8pi sin 01, 
U3 = 8/X2COS02 5 1^4 ~ — 8/i2sin02, 



^0 



17 

-O 



\ 

UJal - iOa2 

Wa2 " Wal / 





-10 0, _ 

I ' ^2 - 






*^3 



/ \ 
1 


\0 -1 0/ 



01 



/O 

1 

Vo 
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-4ri 
z)'' = I -4ri 

4^/2^l -8Ti 




(A4) 



D'^ and Of in the last equation in Eq. (|A3[) are respec- 
tively negative and traceless skew-symmetric matrices. 

With simple calculations, we can obtain the following 
stationary solution of Eq. (|A3[) : 



'm'^{oo) — 0, mf3(cx3) = rn|3(oo) = rn|3(cx3) = 0, 
m^4(oo) = -^pcosidi - 0), m^4(oo) = -^psm{9i - 



f-14 



M = ^{i + ^i-^p2}, 



m,,(oo) = i{l + yr~^^}, 
where p and are given in Eq. ([5]), from which we can 



calculate the stationary state poo and the stationary fi- 
delity F(poo) in Eq. ©. 

Further, recall that the concurrence of the quantum 

state 



/a w 
b z 



z c 
w* d 



can be analytically solved as [52] : 

C(p) = 2 max{|w| - Vb^, \z\ - V^, 0}, 

from which we can obtain the stationary concurrence 
C7(poo) in Eq. dSl). 
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